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ABSTRACT 

We present new X-ray data taken from the ROSAT PSPC pointing archive for 21 
clusters in the Einstein Extended Medium Sensitivity Survey (EMSS). We have sup- 
plemented these data with new optical follow-up observations found in the literature 
and overall, 32 of the original 67 z > 0.14 EMSS clusters now have new information. 
Using this revised sample, we find no systematic difference, as a function of X-ray flux, 
between our measured X-ray cluster fluxes and those in the original EMSS (~ 30% scat- 
ter). However, we do detect a marginal correlation between this observed difference in 
the flux and the redshift of the clusters, with the lower redshift systems having a larger 
scatter by nearly a factor of two. We have also determined the X-ray extent of these 
re-observed EMSS clusters and find 14 of them have significant extents compared to the 
ROSAT Point-Spread Function. Combining these data with extended clusters seen in 
the original EMSS sample, at least 40% of z > 0.14 clusters now have an observed x-ray 
extent thus justifying their classification as X-ray clusters. 

Using our improved EMSS sample, we have re-determined the EMSS X-ray Clus- 
ter Luminosity Function as a function of redshift. We have removed potential mis- 
classifications and included our new measurements of the clusters X-ray luminosities 
and redshifts. We find similar luminosity functions to those originally presented by 
Henry et al. (1992); albeit with two important differences. First, we show that the orig- 
inal low redshift EMSS luminosity function is insufficiently constrained. Secondly, the 
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power law shape of our new determination of the high redshift EMSS luminosity function 
(z = 0.3 — > 0.6) has a shallower slope than that seen by Henry et al. . We have compared 
our new EMSS luminosity functions with those recently derived from nearby sample of 
X-ray clusters and find that the overall degree of observed luminosity function evolution 
is mild at best. This is a result of the shallower slope seen in our EMSS high redshift 
luminosity function and a more robust low redshift determination of the X-ray cluster 
luminosity function from the literature. 

We have quantified the degree of evolution seen in the X-ray cluster luminosity using 
several statistical tests. The most restrictive analysis indicates that our low and high 
redshift EMSS luminosity functions are statistically different at the 95% level. However, 
other tests indicate that these low and high redshift luminosity functions only differ 
by as little as la. These data are therefore, consistent with no evolution in the X-ray 
cluster luminosity function out to z ~ 0.5. 

Subject headings: cosmology:observations - galaxies:clusters:general - galax- 
ies:evolution - surveys - X-rays:galaxies 
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1. Introduction 

One of the most striking results in observational 
cosmology over the past five years has been the re- 
port of rapid evolution in the space density of X- 
ray clusters out to z ~ 0.5. The strongest evidence 
for i such evolution has been presented by Sfoia ex 



al. 



990a ft He nry et al. 19 92 (H92) using a sample 



of clusters serendipitously detected in Einstein IPC 
pointings (the Einstein Extended Medium Sensitiv- 
ity Survey: EMSS). They report that the luminosity 
function ot A-ray selected clusters steepens signifi- 



cantly as a function of redshift with the most distant 
cluster sample having nearly an order of magnitude 
less X-ray luminous clusters (> 10 45 ergs -1 ) than the 
present epoch. 

In recent years, observations by the ROSAT X- 
ray satellite have focused on extending the redshift 
baseline of these evolutionary studies and overall, the 
results qualitatively agreed with the original findings 
of |I|Ienry et al. 1992. For example, all the specific deep 



PSPC pointings towards known optically rich, distant 
clusters of galaxies (z > 0.6) have yet to discover a 



single X-ray bright cluster (Nichol et al. 1994, Cas- 



bander et al. 1994 & Bower et al. 1994). These studies 



however, may be biased since most of their original 
targets were drawn from optically selected catalogues. 

In the long term, such bias can be removed by com- 
piling serendipitous samples of distant X-ray clus- 
ters from the large archive of ROSAT PSPC point- 
ings in the same vein as the EMSS. Several groups 
have already begun this task and their combined ef- 
fort should provide a powerful database for quanti- 
fying the degree of observed X-ray cluster evolution 



(Castander et al. 1995 




Rosati et al. 1995 


, Burke et 


al. 1995, Bcharf et al. 1997 


& Romer et al. 1997 


). This 



work is still in its infancy, mainly because of the large 
amount of optical follow-up needed to verify the pres- 
ence of a distant X-ray cluster and to exclude possi- 
ble contaminants like stars and AGNs. Preliminary 
results from these ROSAT surveys either appear to 



be consistent with the EMSS sample (see Castander 



et al. 1997) 



et al. 1995), or, advocating no evolution (see Collins 



The theoretical consequences of X-ray cluster lu- 
minosity evolution are significant, suggesting that the 
large-scale structure in the universe evolved hierar- 
chically over a relatively short look-back time. How- 
ever, the present results are inconsistent with simple 
scale-invariant hierarchical clustering models, since 



these actually predict an increase in X-ray bright 
clusters with increasing redshift (Kaiser 1986). Sev- 
eral authors have resolved this discrepancy but at 
the expense of breaking scale-invariance and allow- 
ing significant energy input into the gas at high red- 
shift (Kaiser 1991, Evrard & Henry 1991, Cavaliere 



et al. 1993). Alternatively, the problem could be al- 
leviated by steepening the slope of the initial density 
power spectrum on cluster scales, yet this would be 
inconsistent with popula r Cold Dark Matt e r theories 
of struc ture formation ( Davis et al. 1985| , Henry ei 



al. 1992D 



Considering the underlying theoretical importance 
of cluster evolution, we report here a re-examination 
of the EMSS cluster sample since this database is 
the most robust, well-understood catalogue of X- 
ray clusters out to z ~ 0.5 presently available. The 
EMSS, as a whole, has received intense optical scrutiny 
over the past decade and many of its 83 5 X-ray 
sources now have a secure identification (Jq_ 



jioia e\ 



al. 1990b, Stocke et al. 1991). In the next section, we 



describe new data available from the ROSAT point- 
ing archive on nearly half of the z > 0.14 EMSS clus- 
ters used by H92 in their evolutionary study. These 
new data therefore, provide an important check since 
many of these distant EMSS clusters were only ~ 5 
sigma detections in the Einstein data, i.e. there is a 
large uncertainty on their X-ray fluxes. The ROSAT 
re-observations of these clusters represent a signifi- 
cant improvement on this for three reasons: First, the 
clusters received, on average, longer exposure times; 
second, ROSAT has a much lower particle background 
thus improving its sensitivity to low surface bright- 
ness, extended objects; third, the satellite has greater 
on-axis angular resolution. 

This new ROSAT dataset has allowed us to test 
two critical concerns about the EMSS cluster sample: 
i) the accuracy of the X-ray cluster fluxes and there- 
fore, their luminosities; ii) the classification of the 
X-ray emission as a hot, intracluster medium. The 
latter point is very important with regards to iden- 
tication since X-ray point sources such as stars and 
AGNs dominate the X-ray source counts. Therefore, 
only a small error in classifying these objects as clus- 
ters, and vice versa, will make a large difference to the 
final X-ray cluster subsamplc. Therefore, the most 
secure signature of cluster emission is X-ray extent. 
This re-examination of the EMSS is part of a large 
project aimed at searching the ROSAT data archives 
for serendipitous observations of distant X-ray clus- 
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ters; Serendipitous High-redshift Archival ROSAT 



me 



Cl uster (SHAR C) survey (see [Burke et al. 1995| |U1 
et al. 1995j) 



Throughout this paper we assume 
Mpc" 1 and q Q — 0.5 to be consistent 



50kms 
with previous work. 

2. New Data and Analysis 

We present, in Table 1, all the z > 0.14 EMSS 
clusters, from Table 2 of H92, for which new data 
are available including re-observations by the ROSAT 
satellite and/or further optical identification work by 
Gioia & Luppino (1994; GL94) and others. Column 
1 gives the EMSS name of the clusters followed by, 
in column 2, the ROSAT pointing identification num- 
ber. Column 3 presents the measured redshift of the 
clusters which were taken from GL94 or Carlberg et 
al. (1996). Columns 4, 5 & 6 give the results of our X- 
ray source analysis and indicates whether we find the 
cluster X-ray emission extended or not (the number- 
ing given in column 6 is the same as that in Figures 1 
and 2). Column 7 gives the PSPC off-axis angle (in 
pixels) of the clusters and indicates cases where the 
cluster was the original requested target of ROSAT; 
~ 85% of our ROSAT EMSS cluster re-observations 
were again accidental. Column 8 is the averaged 
PSPC exposure time (vignetted corrected), while col- 
umn 9 gives the net source counts (within the band- 
pass 0.4 -> 2.0 keV). Many of the re-observed EMSS 
clusters are detected to significantly higher signal-to- 
noises in our data than that used by H92. Columns 10 
fell are the measured source and background X-ray 
count rates and the final four columns are the X-ray 
fluxes and total luminosities computed from these ob- 
served count rates in both the ROSAT and Einstein 
bandpasses. 

The PSPC pointing data used in this paper were 
obtained from the GSFC ROSAT data archive and 
reduced using a suite of programs we have developed 
as part of the SHARC survey. We refer the reader 



to our forthcoming paper, Romer et al. 1997, that de- 
scribes in detail our data reduction techniques. How- 
ever, for completeness, we briefly highlight the main 
features of our analysis here. The raw X-ray data 
were initially reduced using the methodology outlined 
by Snowden et al. 1994 . For each pointing, an indi- 



vidual exposure map was constructed from the satel- 
lite's orbital information and divided into the binned 
X-ray photon data (15 arcsecond pixels). We concen- 
trated on the "hard" ROSAT band, 0.4 -> 2.0 keV, 



since this substantially reduces the background count 
rate. The exposure and vignetted corrected images 
were then analysed using our own source detection 
algorithm which is based on a mexican-hat wavelet 
transform. This particular wavelet is well-suited to 
analysing PSPC data since it is the second derivative 
of a Gaussian and therefore, this underlying function 
matches closely the (radially averaged) shape of the 
ROSAT Point Spread Function (PSF) which can be 
approximated by a Gaussian. 

A crucial part of our re-analysis of the EMSS sam- 
ple was to determine if the observed X-ray cluster 
emission was extended or not. For this reason, we 
present this particular stage of the analysis in greater 
detail. The extent of all sources in a pointing was 
determined directly from our wavelet data taking ad- 
vantage of the fact that the width of the mexican 
hat filter, at it's zero-crossing points, is equal to the 
width (between inflection points) of the original Gaus- 
sian it was constructed from. Therefore, we grew 
each source outwards from its initial centroid until 
we reached zero in wavelet-space. This source bound- 
ary, in wavelet-space, is then a measurement of the 
size of that source and is simply the product of the 
size of the mexican-hat (we define this CIS (7 yj which 
is the dispersion of the underlying Gaussian) and the 
size of the observed source, which in most cases will 
be the PSF. A final centroid analysis is then carried 
out on each candidate source, using only the pixels 
within this zero-crossing boundary, and their major 
and minor axes determined. 

This procedure was independently carried out for 
three different mexican-hat sizes; a w = 3, 6 and 9 
pixels. The power of this approach is that the differ- 
ent size filters are sensitive to different regimes of the 
off-axis PSF. The smallest filter is best suited to the 
inner region of the PSPC field-of-view, while further 
off-axis, the larger filter sizes become more powerful 
as the PSF increases. For example, at large off-axis 
angles, the smallest filter tends to break-up sources 
because the surface brightness profiles of sources flat- 
tens due to increases in the PSF and reduction in the 
ROSAT telescope effective collecting area. The a w = 
6 and 9 pixel filters help to recover from this. There- 
fore, by changing the filter size in our analysis, we 
naturally accommodate for a changing PSF thus al- 
lowing us to probe further out into the ROSAT field- 
of-view. 

Figures 1 and 2 shows the distribution of the ma- 
jor and minor axes, as a function of off-axis angle, 
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for 1637 high signal-to-noise sources detected in the 
a w = 3 wavelet maps of 407 deep ROSAT PSPC 
pointings. The size of the these sources is the convo- 
lution of the wavelet size with the intrinsic source size 
and therefore, they are larger expected. These point- 
ings have been reduced as part of the SHARC survey 
which is presently investigating all extended sources 
seen in these data. The other a w maps are very simi- 
lar and are therefore, not shown. The only difference, 
as highlighted above, is the area over which the PSF 
locus is well-defined has been shifted to higher off-axis 
angles. Clearly visible in Figure 1 is the expected in- 
crease in the ROSAT PSF with off-axis angle above 
which we must determine if a source is extended or 
not. Also shown is a line which represents sources that 
are > 3er extended and was empricially determined, 
as a function of off-axis angle, from these sources (see 
Ftomer et al. 1997 for further information). 

The off-axis angle and major and minor axes of all 
EMSS clusters detected in the 3 pixel wavelet map 
are presented in Table 1 (column 6) and are plotted 
in Figures 1 and 2. An EMSS cluster is designated 
as extended if its major and/or minor axis is above 
this 3<7 line in any of the 3 wavelet maps, for example, 
MS0811. 6+6301 was selected as extended based solely 
on the a w = 6 wavelet information. A question mark 
in Table 1 signifies clusters that do not satisfy these 
criteria, but appear to lie away from the point-source 
locus in two or more of its wavelet representations. 

The X-ray count rates of the EMSS clusters were 
measured interactively using the IRAF/PROS soft- 
ware package. This differs from the approach used 
in the SHARC survey and therefore, we describe it 
fully here. In the original EMSS, the cluster count 
rates were computed within a fixed angular aperture 
(2.4 x 2.4 arcmins) and were later corrected, via a 
cluster model, to the total X-ray flux. We, however, 
had the freedom to select a fixed metric aperture thus 
scaling it sppropriately with redshift and the PSF for 
each individual cluster (we convolved the angular size 
of the metric aperture with the expected size of the 
PSF). For aa assumed King profile with a core radius 
of 0.25Mpc and f3 = |, we computed the angular size 
of an aperture for each cluster that would encompass 
85% of the total cluster flux, irrespective of redshift 
and off-axis angle (the fluxes and luminosities given 
in Table 1 are corrected for the 15% of light expected 
outside the aperture). Furthermore, the interactive 
approach allowed us examine the environment of each 
cluster and exclude any nearby sources (both in the 



background and source apertures). Wherever possi- 
ble, the background was estimated from an annulus 
around the source aperture having a similar number of 
pixels as the source aperture. If this was not possible, 
a nearby area devoid of sources was used, ensuring 
that it was at the same off-axis distance as the source 
so the vignetting corrections would be comparable. 

These measured count rates were finally converted 
to fluxes (erg cm~ 2 s _1 ) by integrating a 6keV Raymond- 
Smith (RS) spectrum through the ROSAT response 
function and accounting for absorption using HI col- 
umn densities interpolated from Stark et al. (1992) 
and the model by Morrison & McCammon (1983). 
This therefore, gave us the X-ray cluster flux in the 
ROSAT bandpass outside our own Galaxy. We then 
converted the ROSAT fluxes to that expected in the 
Einstein bandpass again using a 6keV RS spectrum. 
We chose this model cluster spectrum to be fully con- 
sistent with the H92 analysis. The cluster luminosi- 
ties (including k-corrections) for both the ROSAT 
and EMSS bandpasses are given in columns 14 & 15 
of Table 1 respectively. 

Ideally, we would like to convert our ROSAT count 
rates to Einstein fluxes using the individual spec- 
tra of these EMSS clusters. However, for a vast 
majority of the EMSS such data does not exist. 
Therefore, to quantify this potential systematic effect, 
we re-computed the expected Einstein fluxes of all 
our ROSAT EMSS detections using the Luminosity- 



Temperature relationship of Wang & Stocke 1993 
and iterated it until we gained a stable solution 
(within 0.1% between iterations). This usually took 
< 4 passes and always changed the final iterated 
cluster luminosity by less than 6%. For example, 
MS0015.9+1609 shifted from L x (44) = 22.78 ergs" 1 
- which is given in Table 1 - to L x (44) = 23.72 ergs -1 
for an iterated temperature of 11.35keV. Our count 
rate to flux conversion was also relatively insensitive 
to variations in the RS metal abundance and varied 
by less than 10% over a reasonable range of metallic- 
ities. 

In Figures 3 and 4, we present all available ROSAT 
data on the distant z > 0.14 EMSS clusters dis- 
cussed in this paper. The size of the individual plot- 
ting windows in Figure 3 is equal to the diameter of 
the aperture used to compute the cluster count rates. 
Also shown is the measured position angle of the ob- 
served x-ray emission and the radial direction towards 
the center of the PSPC. In Figure 4, we present the 
ROSAT radial X-ray surface brightness profiles for 
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the clusters as well as the appropriate off-axis ROSAT 
PSF. For the PSPC data, this was computed from the 
point-source locii seen in Figures 1 and 2 (not the 30- 
line) . For the HPJ data, we used the model of David 
et al. (1993) for the PSF. These plots clearly show 
that many of the re-observed EMSS clusters are high 
signal-to-noise X-ray detections and are extended. 

3. Comparison of Cluster X ray Flux Mea- 
surements 

One of the motivations of our work was to check 
the robustness of the EMSS flux measurements. In 
H92, the total flux was computed, using a King model, 
from the detected flux within a fixed angular aper- 
ture. In contrast, we attempt to measure the total 
flux directly using a fixed metric aperture. For the 
21 ROSAT-EMSS clusters, Figure 5 presents the dif- 
ference between our ROSAT mesurement of the to- 
tal flux of the EMSS clusters in the Einstein band- 
pass (column 13 in Table 1) and that given by H92, 
versus our total ROSAT-EMSS flux. The error bars 
plotted in Figure 5 were derived from the propaga- 
tion of the errors published on the net counts by H92 
and in Table 1 of this paper (column 9) . Using a x 2 
test, we have compared the data in Figure 5 with the 
expected relationship i.e. the zero line in Figure 5. 
For all the data points, we find x 2 — 57 for 20 de- 
grees of freedom (d.o.f) which is highly improbable. 
However, much of this disagreement is due to the two 
points at the highest total flux (MS0451. 5+0250 and 
MS0735. 6+7421). Therefore, if we constrain our anal- 
ysis to F x < 2 x 10 -12 ergs -1 cm 2 and F x < 10~ 12 ergs" 
we find x 2 = 33 (18 d.o.f) and X 2 = 12 (12 d.o.f) re- 
spectively. These indicate that the observed scatter 
(28%) at lower flux levels is consistent with the ob- 
served errors on the two datasets (ours and H92). 

A more interesting plot to make is to plot the flux 
differences in Figure 5 as a function of redshift. This 
is shown in Figure 6. In this case, there does appear 
to be a noticeable decrease in the flux differences with 
increased redshift. This is quantified by splitting the 
sample at z = 0.3, with the lower redshift clusters 
having a mean absolute difference of 39%, while for 
the higher redshift data, the mean absolute difference 
is only 19%. The standard error on both these means 
is ~ 30%, since both were determined using ~ 10 
clusters. Therefore, this observed drop in the mean 
absolute flux difference is marginally significant. Such 
a decrease with redshift is consistent with the H92 
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Fig. 1. — The major and minor axes for 1637 high 
signal-to-noise sources detected in 407 deep ROSAT 
PSPC pointings versus their off-axis angle (the x-axis 
the convolution of the source size with the wavelet 
filter size). The change in the PSPC PSF with off- 
axis angle is visible and was defined empiricially (the 
dotted lines represent > 3<r extended sources). The 
star symbols are z > 0.3 EMSS clusters and have the 
same numbering as given in column 4 of Table 1. The 
unfilled circles are the 2 EMSS clusters that only have 
lower limits measured for their extents. These sources 
were found in the 3 pixel wavelet map. 
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Table 2: Maximum likelihood Parameterisation of XCLFs derived from the original EMSS Data Sample 
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methodology since at higher redshifts there would be 
a greater amount of the clusters' total flux within the 
fixed angular aperture compared to lower redshifts, 
thus the required correction would be smaller. 

4. Measurement of Extended Cluster X-ray 
Emission 

It is imperative that the classification of all EMSS 
clusters should be made as secure as possible. This 
is straightforward using the measured extent of the 
observed X-ray emission from the cluster. Figures 
1 and 2 present the observed extent of our ROSAT- 
EMSS clusters against that expected for X-ray point- 
sources (as a function of off-axis angle). Overall, 14 
of the 21 clusters with new ROSAT data are flagged 
extended in Table 1 (column 4). For 3 of these ex- 
tended sources, H92 also detected an extent in the 
Einstein data (these clusters are highlighted in Table 
1). By combining our extent results with H92, a to- 
tal of 25 z > 0.14 EMSS clusters now have sufficient 
X-ray information to measure a significant X-ray ex- 
tent (40% of the whole EMSS sample). Constraining 
ourselves to z > 0.3, this figure drops to 30%. For 
these extended clusters, therefore, we are most prob- 
ably seeing X-ray emission from a hot, intracluster 
gas which undoubably strengthens their classification 
as X-ray clusters. 

Next we discuss the clusters marked as non-extended, 
or uncertain, in our database. Some of these are the 
result of low signal-to-noise, low intrinsic X-ray lu- 
minosity and/or uncertainties in the PSF at large off- 
axis angles. For example, MS0418. 3-3844 is very close 
to the ROSAT rib structure resulting in an unrealis- 
tic major and minor axis in all 3 wavelet maps and 
therefore, we only present a lower limit to the clus- 
ter size in Table 1. The same problem occurred for 
MS81 1.6+6301 in the a w = 3 wavelet map, but it was 
observed to be extended in the larger wavelet maps. 

For EMSS clusters MS1209.0+3917, MS1219. 9+7542, 
MS1512.4+3647 and MS2137.3-2353, there exists sig- 



nificant archival ROSAT HRI data which can help 
classify these systems. In the PSPC data, they all 
lie on - or near - the observed point source locus 
in all three wavelet maps (see Figures 1 and 2). 
They also appear unextended in Figure 4 Consider- 
ing the length of exposure time gained in these PSPC 
data, we suggest that these clusters either have com- 
pact cores, or, that they are mis-identifications. For 
MS1209.0+3917 and MS1219. 9+7542, the HRI indi- 
cates that these X-ray sources are point -like with lit- 
tle evidence for any X-ray extent. The optical CCD 
image of MS1209.0+3917by GL94, and their own sus- 
picions, supports the idea that this cluster is a mis- 
identification. For MS1219. 9+7542, this system is a 
poor group with a dominant galaxy showing mod- 
erate [Oil] emission in its' optical spectrum (GL94 
and Stocke et al. 1991). The X-ray and optical data 
combined, appears to indicate that the x-ray emission 
for these two EMSS clusters may be due to an AGN 
and not a hot, intracluster medium. We therefore, 
remove these clusters in any further analysis of the 
EMSS sample. 

For MS2137.3-2353 and MS1512. 5+3647, the HRI 
data indicates that the X-ray emission from these 
systems is extended compared to the PSF. We fit- 
ted a King model to the HRI data and measured 
core radii of 17 ± 8 (j3 = 0.77 ± 0.28) and 16 ± 10 
(/3 = 0.63 ± 0.23) arcseconds for MS2137.3-2353 and 
MS1512. 5+3647 respectively. Clearly, both clusters 
could not have been expected to be resolved by the 
PSPC instrument since even on-axis, the PSPC PSF 
has a ~ 20 arcseconds FWHM (illustrated in Fig- 
ure 4). For both systems, the angular extent of the 
X-ray emission represents a metric core radius of less 
than 100 kpc; significantly smaller than the rest of the 
EMSS sample. For MS2137.3-2353, the lensing work 
of Mellier et al. (1993) agrees with this metric core 
radius, while the optical appearance of this system 
is dominated by a single, large cD-like galaxy with 
a similar metric size as the observed X-ray emission 
(GL94). For MS1512. 5+3647, Stocke et al. (1991) has 
classified this system as a "cooling-flow galaxy" which 
is defined to be a galaxy with similar optical char- 
acteristics as the dominant galaxy in a cooling flow 
clusters but which is no associated with a rich, opti- 
cal cluster of galaxies (see the CCD image in GL94). 
Carlberg et al. 1996 has extensive redshift informa- 
tion on this EMSS cluster and measure a velocity dis- 
persion of a v = 690 km s -1 . However, this system 
appears to have substantial redshift contamination 



(as seen in Figure 1 of their paper) and may indicate 
that MS1512. 5+3647 is a projection effect. More data 
are required on these two systems to unambiguously 
classify them. We do not remove MS2137. 3-2353 and 
MS1512. 5+3647 in our later analysis, but note that 
our results were quite insensitive to their inclusion, 
or, removal. 

No high quality ROSAT HRI data exists for the re- 
maining 2 clusters marked as non-extended, or uncer- 
tain, in Table 1 (MS12O8.7+39280& MS1617.1+3237). 
For these two, we therefore performed Monte Carlo 
simulations to assess whether we had the required 
signal-to-noise to observe an extent or not (these 
simulations are presently being developed as part of 
the SHARC survey to help characterise the selec- 
tion function of that survey and will be presented 
in detail in a forthcoming paper). Briefly, the sim- 
ulations used here involved adding false clusters to 
the relevent ROSAT PSPC pointings (Table 1) us- 
ing an underlying King profile of core radius 0.25 
Mpc and /3 = |. These false clusters were given the 
same signal-to-noise as the real real detections and 
were positioned at the same off-axis angle (usually 
opposite the real cluster in the PSPC field-of-view). 
For each cluster, 10 iterations were constructed and 
passed through our source detection software. This 
allowed us to determine the frequency with which 
these false sources were flagged as extented. Using 
these simulations as a guideline, it would appear that 
MS1208. 7+3928 is more compact than expected since 
it was always detected as extended in the simula- 
tions. MS1617. 1+3237, however, remains uncertain 
since the simulations indicated that we do not have 
high enough signal-to-noise to unambigously detect 
an extent given this nominal profile (it was flagged as 
extended only ~ 70% of the time) . 

5. The X ray Cluster Luminosity Function 
5.1. A Re Determination of the H92 Result 

As mentioned earlier, the most striking result to 
come out of the EMSS cluster sample is evidence for 
rapid negative evolution in the X-ray Cluster Lu- 
minosity Function (XCLF; H92). Considering the 
amount of new information gained on the EMSS sam- 
ple of clusters over the past few years, we have revis- 
ited the H92 result here. Our motivation for this are 



2 MS1208.7+3928 was not detected in the ROSAT HRI pointing 
701844 towards MS1209.0+3917. 
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Fig. 3. — A greyscale plot of the count rates of the 21 EMSS cluster with new ROSAT data given in Table 1. 
North is upwards with East to the left. The solid line represents the measured position angle of the cluster and 
is only computed for clusters with minor/major > 0.8, thus avoiding circular clusters. The numerical value - in 
degrees - of the position angle is the top number shown in the top-left hand corner of the window. The dashed 
line is the direction towards the center of the PSPC and helps the reader determine if the asymmetrical PSPC 
PSF has influenced the determination of the position angle. In some cases, this looks likely The dashed line is not 
plotted if the cluster was nearly on-axis. The count rate data are scaled linearly between zero and the maximum 
observed count rate (these maps have been vignetted and exposure corrected). The contribution of the PSF has 
not been removed and these images have also been lightly smoothed with a Gaussian of width one pixel to reduce 
the binning noise. The size of the plotting window is given as the bottom number in the top left-hand corner and 
is equal to the aperture used in the flux determination. Other sources near the cluster were removed prior to any 
flux determination i.e. MS1208. 7+3928. 
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two-fold: First, to determine if actual changes in spe- 
cific clusters make a significant systematic difference 
to the observed XCLF; secondly, to assess the robust- 
ness of the result to small, yet significant, changes 
in the overall cluster population. For instance, if the 
XCLF was radically altered, in any sense, by changes 
to a few clusters, it would indicate that the sam- 
ple was not stable and remove much of its predic- 
tive power. In conjunction, we have re-assessed the 
statistical significance of X-ray evolution originally 
reported by H92. 

We initially re-computed the XCLF using the same 



1/V a test as employed by H92 (Avni fc Bahcall 1980) 



This involved computing the maximum observable 
volume (V a ) for each cluster, 



N 



max ; ^rain ) j 



(1) 



where ili is the area of the EMSS survey to a given 
sensitivity limit (see Table 3 in H92), N is the total 
number of different sensitivity limits (19 in total from 
Table 3 in H92) and z m i n is the lower redshift limit 
for a given redshift shell. For z max , the smaller of 
either the maximum detectable redshift for a given 
cluster, or, the high redshift limit in a given redshift 
shell was used. The maximum detectable redshift was 
computed in the same fashion as in H92 (Eqn. 2 in 
their paper) using the sensitivity limits in Table 3 
of H92 and either our new ROSAT determination of 
the cluster luminosity (using the same csomology and 
k-corrections) or that given in Table 2 of H92 (for 
clusters with no new ROSAT information). 

These individual cluster volumes were then con- 
verted into a XCLF by binning them as a function of 
luminosity according to 



1 



V a .j AL' 



(2) 



where n is the number of clusters in the luminosity 
bin AL wide. Again, we use new ROSAT measure- 
ments of the clusters X-ray luminosities where appro- 
priate. Before we discuss our new determinations of 
the EMSS XCLFs, we note that we reproduced the 
binned H92 XCLFs as shown in Figure 2 of their pa- 
per using only the data given in Table 2 of H92. The 
re-determinations of the H92 XCLFs are shown in 
Figure 7 for comparison and serve as a good check of 
our methodology. 



Meanwhile, Figure 8 shows our new determinations 
of the EMSS XCLFs in the same three redshift shells 
as used by H92; z=0.14 to 0.2, 0.2 to 0.3 and 0.3 to 
0.6. We have used, where appropriate, the new data 
presented in Table 1 of this paper, which includes 
removing clusters that we have flagged as either un- 
extended or uncertain (see the discussion above and 
Table 1 where we have marked removed clusters with 
a *). For clusters with no new information associated 
with them, we simply used the original fluxes and 
luminosities quoted in H92. Poisson error bars are 
shown and were obtained from the confidence limits 



tabulated by |Gehrels 1986 . 

Clearly, the three new XCLFs presented in Fig- 
ure 8 are similar to those shown in H92 and Fig- 
ure 7. This indicates that the EMSS sample is in- 
ternally robust, since 50% of the clusters used in our 
new EMSS XCLFs have either a changed classification 
(cluster or not), a new flux measurement (different 
by up to 60%), and/or a new redshift. Certain clus- 
ters have changed substantially in luminosity, popu- 
lating different bins in our XCLF than they did in the 
H92 result. This is certainly true in the highest red- 
shift shell where the 2 clusters MS0015. 9+1609 and 
MS0451. 6-0305 have been elevated in total luminos- 
ity to occupy a previously unoccupied bin centered 
at 2.5 x 10 45 ergs _1 . We note here that other re- 
determinations of the ROSAT flux for these two clus- 
ters are in good agreement with ours (Hughes, private 
communication, Neumann & Bohinger 1996, Donahue 
& Stocke 1995). 

However, these binned XCLFs do not allow us 
to judge the statistical significance of our work and 
therefore, we parameterised our luminosity functions 
using a power-law luminosity function of the form, 
71(^44) = K £44° , where L44 is the clusters luminos- 
ity in units of 10 44 ergs -1 (H92). For every cluster in 
a particular redshift shell, we compute the probabil- 
ity that we would see this cluster, at that luminosity, 
given the sensitivity limits of the EMSS (in H92 Ta- 
ble 3) and the model. This probability (P) can be 
expressed in full as, 



P(L Z ) = 1 f V(L) L a 6(1., - L) dL, 



(3) 



where Li is the observed luminosity of a cluster (taken 
from Table 1 of this paper or Table 2 from H92), V(L) 
is the available volume within which that cluster could 
have been seen and 8{Li — L) is a delta function which 
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Fig. 4. — The radially-averaged profiles for 21 EMSS clusters with new ROSAT PSPC data (given in Table 1) and 
4 EMSS clusters with available HPJ data. The x-axis is in arcseconds from the intensity-weighted centroids of the 
clusters, while the y-axis is the observed X-ray surface brightness (counts s -1 arcsec~ 2 ). The binsize used for the 
PSPC and HRI data was 20 and 6 arcseconds respectively. The dashed lines show the expected size of the PSF at 
that given off-axis angle (as given in Table 1). The PSPC PSF was empirically determined from the data given in 
Figures 1 and 2 (we use the observed point-source locus not the > 3tr line shown). The HRI PSF was taken from 
David et al. 1993 
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Table 3: Maximum likelihood Paramctcrisation of XCLFs derived from our New EMSS Data Sample 



Rcdshift shell 


Number 


a 


K(10 7 Mpc-' 6 L^ L ) 


0.14- 


0.2 


19 


2.44 ±0.34 (2.19 ±0.21) 


7.9 ±1.8 (5.85 ±0.25) 


0.2 - 


0.3 


21 


2.66 ±0.30 (2.67 ±0.26) 


6.5 ± 1.7 (6.82 ±0.51) 


0.3 - 


0.6 


21 


3.03 ±0.30 (3.27 ±0.29) 


7.9 ±3.0 (12.33 ±3.87) 


0.14- 


0.3 


40 


2.62 ±0.22 


7.2 ± 1.3 


With Luminosity Errors 


0.14- 


0.2 


19 


2.60 ±0.37 


8.7 ±2.0 


0.2 - 


0.3 


21 


2.73 ±0.31 


7.0 ± 1.8 


0.3 - 


0.6 


21 


3.08 ±0.32 


8.5 ±3.4 



T 




5 10 15 



Y i i . i i i . i i 

20 40 60 

ROSAT Total Flux (10" 13 erg/s/cm 2 ) 

Fig. 5. — The percentile difference between our es- 
timate of the EMSS total flux (derived from our 
ROSAT data) and the total EMSS flux quoted by 
H92, as a function of our estimate of the EMSS total 
flux. The errors plotted were derived from the com- 
binations of the errors on the net counts given in H92 
and our errors in Table 1. 




Redshift 

Fig. 6. — The same as Figure 4, but plotted as a 
function of EMSS cluster rcdshift. 
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is unity at Li and zero for all other values of L. This, 
therefore, assumes no error on the observed luminosi- 
ties (see below) and Eqn. 3 can be reduced to 

P{L i ) = ^V a L^AL, (4) 

where V a is given in Eqn. 1 (the available search vol- 
ume of the cluster) and AL is chosen to be vanishingly 
small so only a L t cluster can fall within this interval 
(see Cash 1979 for full details). In Eqn. 3 and 4, N 
is the normalisation such that 

/>oo 

N= V(L)L- a dL (5) 

intergrated over the entire same range of luminosities. 
The likelihood of a given model is therefore, 

n bj 

C(a) = [] P(Li\a) (6) 

using Equations 3,4 and 5. 

Since our ML method is slightly different from the 
approach implemented by H92, we ensured that our 
ML fits to the original H92 data were compatible with 
those published in their paper. In other words, we at- 
tempted to fully reproduce the H92 result using only 
the data in Table 2 of their paper. We summarise 
these results in Table 2 and the fits are shown in Fig- 
ure 7. Overall, we obtained excellent agreement with 
the H92 results with only one significant discrepancy; 
we see a 50% large error on the low redshift XCLF 
than H92 (In Tables 2 and 3, we present in paren- 
theses the slopes, and errors, on all 3 EMSS XCLFs 
as published in H92). The source of this discrepancy 
remains unclear but may be the result of us have a 
different likelihood distribution, for this low redshift 
shell, than H92 and Gioia et al. (1990a). This would 
effect our error analysis because we integrated the 
observed likelihood distribution directly, while H92 
and Gioia et al. (1990a) used an analytical approach 
which assumes the likelihood distribution is a Gaus- 
sian (Gioia, private communication). We performed 
this latter approach on our data and did find a smaller 
error on the slope of the low redshift EMSS XCLF 
compared to that presented in Table 2. However, it 
was not large enough to fully explain this discrepancy 
between us and H92. 

Table 3 summaries the ML parameterisations of 
our new EMSS XCLFs and includes (in parentheses) 





CO 

1 


— i — r 


i 1 i i i i 1 i i i 


i i i 1 1 


- 




CO 

1 








- 




o 

T — 1 

1 


_0.14 < z < 0.2 

i i i i 1 i i i i 1 i i i 


i I i 1 1 




1 


CO 




i ii i i i ii i i 




+- 


< 


1 


< 








CO 
1 












(Mpc" 


CO 

1 










log N 


o 

1 
1 


_ 0.2 < z < 0.3 
■ i i i l i i i i l i i i 


\. < 
i I i 1 1 


> 




CO 

1 




i ii i i i ii i i 




+- 




CO 

1 












o 

1 


_ 0.3 < z < 0.6 
i i i i l i i i i l i i i 


i i i 1 1 


i 



43.5 44 44.5 45 45.5 
log 1^ (0.3 - 3.5 keV) 



Fig. 7. — Binned X-ray Cluster Luminosity Functions 
in three redshift shells. These were computed by us 
using only the original data in H92 (we have used 
different bin centers to those used by H92). The solid 
line represents our best ML fit while the dotted line 
is the best ML fit from H92. 



13 




J I I I I I I I I I I I I I I I I I L 



43.5 44 44.5 45 45.5 
log 1^ (0.3 - 3.5 keV) 

Fig. 8. — Binned X-ray Cluster Luminosity Functions 
in three redshift shells. In this case, we have used the 
improvements detailed in Table 1 of this paper in con- 
junction with data from the original H92 data. The 
solid line represents our best ML fit as seen in Table 2. 
The square symbol in the bottom plot represents the 
volume density of very high redshift EMSS clusters 
discussed in the text. 



the original H92 results for comparison. The normal- 
isation of the power-law XCLF (K) was computed 
by requiring Eqn. 5 to give the observed number of 
clusters. Our ML parametric representations are in 
reasonable agreement with those published by H92. 
This should not be too surprising given the visual 
similarity in the binned XCLFs. The largest discrep- 
ancy is in the slope of the highest redshift XCLF, 
where we find a shallower slope than H92 by ~ la. 
This change is almost exclusively due to the increased 
dynamic range in luminosity we are probing since we 
have a new highest luminosity bin compared to H92. 
We also find that our determination of the low red- 
shift XCLF has a steeper slope compared to the orig- 
inal H92 result (we were able to re-produce the H92 
slope when we used only their data). This is some- 
what surprising since we have only changed 4 clusters 
in this redshift shell (2 changed in luminosity and 
2 clusters were removed; see Table 1). We have in- 
vestigated this discrepancy by re-computing the low 
redshift XCLF removing each of these 4 clusters one 
at a time. The most significant effect is noticed when 
MS2318. 7-2328 is removed because of a lower redshift 
measurement by Romer (1994). This cluster is the 
highest luminosity cluster in the original H92 deter- 
mination of the low redshift XCLF and therefore, it 
is understandable that its' removal has the largest ef- 
fect on our XCLF re-determination (also it explains 
why we see a steeper slope since we have removed the 
highest luminosity cluster in that shell). Overall, this 
low redshift shell is poorly determined since it is very 
sensitive to small changes in the dataset. 

To quantify the statistical significance of our ML 
fits, we present in Tables 2 and 3 the computed er- 
rors on our ML parameters. These were computed by 
re-normalising the observed likelihood distributions 
(by integrating over the whole distribution and divid- 
ing it by the total) for the three redshift shells and 
determining the interval that enclosed 68% of the val- 
ues. Figure 9 presents our observed likelihood distri- 
butions for the three redshift shells and can be ap- 
proximated by a Gaussian. This figure is extremely 
enlightening, since it shows that the slope (a) of the 
low redshift parametric XCLF is insufficiently con- 
strained since it has a wide likelihood distribution. 
Overall, the high redshift XCLF parameterisation is 
only mildly inconsistent with the low redshift data. 
The probability of not getting the high redshift slope 
- from the low redshift data - is only 93%. This is 
caused by a combination of the low dynamic range 
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- in luminosity - that the low redshift data extends 
over and the changes in the XCLF slopes mentioned 
above. Moreover, the difference in slope between the 
low and middle redshift XCLFs is clearly insignificant 
with the probability of obtaining the middle redshift 
shell maximum likelihood slope - from the low red- 
shift data - being 51%. 

Since the error on the XCLFs is crucial to any mea- 
surement of XCLF evolution, we tested our ML er- 
ror analysis by performed Monte Carlo simulations 
(> 2000 iterations per XCLF). This involved gener- 
ating fake cluster databases - with the same total 
number of clusters as the real data - drawn at random 
from the parameterised H92 XCLFs. For all three red- 
shift shells, the mean of the subsequent distribution 
of fitted XCLF slopes measured for these fake cluster 
catalogues was identical to the original input values. 
Moreover, the distributions were Gaussian-like and 
had dispersions very similar to those presented in Ta- 
bles 2 and 3. Therefore, our error estimates on the 
real XCLFs given in Tables 2 and 3 are consistent with 
the scatter expected for these number of clusters (in 
each redshift shell). 

5.2. Further Tests of XCLF Evolution 

In this section we move beyond the analysis used 
by H92 to incorporate different statistical tests and 
independent XCLF results published over the last few 
years. We also investigated splitting the sample into 
different redshift shells. 

Our analysis has highlighted that the low redshift 
EMSS XCLF is insufficiently constrained. Recently, 
however, Ebeling et al. (1995, 1996) has published a 
new low redshift XCLF using the ROSAT All-Sky 
Survey (RASS) Brightest Cluster Sample comprising 
of 173 z < 0.3 X-ray clusters. From their analy- 
sis, they found that the data was consistent with no 
XCLF evolution out to z = 0.3. This can be seen in 
Figure 10, where we compare this new RASS XCLF 
with our re-determiations of the EMSS XCLFs. We 
have now combined the lower two redshift EMSS 
shells (i.e. z = 0.14 to 0.3) since our analysis indicates 
that there is little statistical difference between the 
two. Furthermore, it helps in our search for evidence 
for XCLF evolution because it increases the dynamic 
range - in cluster luminosities - studied thus provid- 
ing a greater constraint on the fitted ML parameters 
(see Figure 9 and Table 3). 

Clearly, there is excellent agreement between the 




Fig. 9. — The likelihood distributions for the three 
different redshift EMSS X-ray Cluster Luminosity 
Functions. The dashed line is the combined luminos- 
ity function (z = 0.14 to 0.3) discussed in the text. 
The solid points -with errors- represent the maximum 
likelihood slopes as tabulated in Table 2. 
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44 44.5 45 45.5 

log 1^ (0.3 - 3.5 keV) 

Fig. 10. — Comparison of the low redshift XCLF (solid line) with the high redshift XCLF (dotted line) and the 
XCLF of Ebcling et al. (1995, 1996) computed from the RASS (dashed line). The latter has been shifted to the 
same energy bandpass. The solid symbols (•) are binned high redshift luminosity function, while the open circles 
(o) are the low redshift case. The solid square is the volume density of very high redshift (z > 0.6) clusters and 
the arrow (|) marks the lower limit of Nichol et al. (1994) in this redshift range. 
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Fig. 11. — V/V m ax versus L x for the EMSS sample. 
The left-hand panel is for the combined low (0.14 to 
0.3) redshift shell (•) against the high rcdshift EMSS 
shell (open symbols) . The right-hand panel is for the 
original low (0.14 to 0.2) redshift shell (*) against the 
same high redshift shell data. 



RASS XCLF and our combined lower redshift EMSS 
XCLF. This is very impressive since the two samples 
cover different redshift regimes i.e. the RASS sam- 
ple only has 20% (8%) of its clusters at z > 0.14 
(z > 0.2) - the lower limit of the EMSS data - and 
all having L x > 5 x 10 44 ergs -1 ). The fact that these 
two XCLFs, from different surveys over different red- 
shift ranges, agree so well justifies our original mo- 
tivation for combining the two lower EMSS samples. 
Clearly the emphasis should now be on determining 
the degree of evolution above z = 0.3. We say this ir- 
respective of our own work, since Ebeling et al. (1995, 
1996) sees no evidence for evolution at rcdshifts lower 
than this using a much larger sample of clusters. We 
have therefore, compared this new combined lower 
redshift EMSS XCLF with that derived at high rcd- 
shift. Although the fitted ML slopes of these two 
EMSS XCLFs are different (z = 0.14 to 0.3 compared 
to z = 0.3 to 0.6), a close inspection of Figure 9 and 
Table 3 highlights that this discrepancy is only mildly 
significant. 

To quantify this significance, we have implemented 
a 2D Kolmogorov-Smirnov (KS) test, in the L x - 
V/Vmax plane, comparing the two redshift shells, 0.14 
to 0.3 and 0.3 to 0.6, from our improved EMSS sam- 
ple. The V/Vmax was calculated using the approach 
outlines in Avni & Bahcall (1980) using a uniform 
density distribution. This statistical analysis removes 
the need for any binning, or fitting, of the data and 
directly test the hypothesis that these two distribu- 
tions were drawn from the same underlying parent 
distribution. Therefore, for this application, it is a 
very powerful method. Figure 11 shows the data 
displayed in the L x -V/V max plane for the two red- 
shift shells mentioned above. A KS test of these two 
gives a probability of 15% that they were drawn from 
the same underlying distribution. This drops only 
slightly to 11% if we limit ourselves to comparing 
the z = 0.14 to 0.2 shell to the high redshift shell. 
In these KS tests, we have constrained the data to 
L x (0.3 — 3.5) > 10 44 ergs -1 since below this limit the 
high redshift shell starts to become seriously incom- 
plete (the EMSS does not have the sensitivity to probe 
such clusters at these high redshifts). If we compare 
all the data without a luminosity cut, the above KS 
probabilities - that the two distributions are the same 
- are 9% and 5% respectively. 

For comparison, we also computed the mean V/V max 
in the high redshift shell (over the same range of L x ) 
for both our improved EMSS sample and the original 
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H92 sample. We found 0.44 ± 0.08 for our data and 
0.36 ± 0.06 for H92 (the standard error on the mean 
is quoted). This again indicates that the evidence 
for evolution in the high redshift shell has decreased 
since our measured mean V/V ma x is now consistent 



with the expected value (0.5; Avni & Bahcall 1980). 

Most statistical analyses of XCLFs do not include 
observed luminosity errors. As can be seen in Figures 
5 and 6, this can be a large effect and a function of 
redshift. For completeness therefore, we re-computed 
our ML parameterisations of the XCLFs but replacing 
the delta function in Eqn. 3 by a Gaussian whose 
width is equal to the observed luminosity error (using 
either the error on the flux given in Table 2 of H92 
or our error on the ROSAT flux in Table 1 of this 
paper). This is computationally intensive since for 
each cluster the probability of it being observed is 
now an integral over a range of possible luminosities 
(because of the error in the luminosity) . The search 
volume for each possible luminosity is re-computed 
(Eqn. 1) since z max changes with luminosity. 

In Table 3, we present these new ML parameter- 
isation of the low, middle and high redshift EMSS 
XCLFs. As can be seen, their are two effects in adding 
such errors. First, all the slopes have steepened with 
respect to their previous determinations (with no er- 
rors). Second, the degree to which they have steep- 
ened is different for each redshift shell, with the lower 
redshift shells changing the most. Therefore, the sta- 
tistical significance of any steepening with redshift of 
these slopes is now lessened; the low and high redshift 
shells are only inconsistent at the la level. 

Finally, we turn our attention to the form of the 
very high redshift XCLF i.e. z = 0.6 to 0.9. In the 
past few years, several groups have used the ROSAT 
satellite to extend the study of XCLF evolution into 



this redshift regime (Nichol et al. 1994, Castander 
et il. 1994). Moreover, further optical follow-up by 



GL94 has revealed two new EMSS clusters that have 
z > 0.6 (see also Luppino & Gioia 1995). Although 
the data are still very limited, we have plotted the 
volume density of these two very high redshift EMSS 
clusters as a single point in Figure 10. Also plotted 
is the lower limit for lower luminosity cluste rs in this 
redshift shell as derived by Nichol et al. 1994 . In both 
case, these data are consistent with the lower redshift 
EMSS XCLFs. 

Using the EMSS sensitivity limits given by H92, 
we have computed that the lowest detectable clus- 
ter luminosity in this very high redshift shell is 4.1 x 



10 44 ergss _1 . Given the low redshift XCLF parame- 
terisation (redshift 0.14 to 0.3; Table 3), we estimate 
that the EMSS should have detected a total of 10 
clusters brighter than this limit at these very high 
redshifts (we would expect three of these clusters to 
have L x > 10 45 ). This discrepancy is either evidence 
for evolution in the bright end of the XCLF at these 
high redshifts, or, a reflection of the EMSS classifica- 
tion procedure (the highest redshift clusters are the 
hardest to find and measure). 

6. Conclusions 

We present in this paper a thorough re-examination 
of the EMSS cluster sample using new data from the 
ROSAT PSPC pointing archive. Furthermore, we in- 
clude the latest information on the optical follow-up 
of this sample of clusters. In total, 32 of the origi- 
nal 67 z > 0.14 EMSS clusters used by H92 to study 
X-ray cluster evolution have new information associ- 
ated with them, of which, 21 have new X-ray data 
obtained from the ROSAT archive. For these clus- 
ters, we find no systematic difference, as a function 
of X-ray flux, between the original EMSS flux esti- 
mates and our ROSAT re-determinations (28% scat- 
ter). Our analysis did however, indicate that this 
scatter may be correlated with cluster redshift with 
the high redshift clusters having a lower scatter by a 
factor of two. This is consistent with the expected 
trend given that the original EMSS fluxes were com- 
puted within a fixed angular aperture. 

We determined the extent for 21 of the ROSAT 
re-observed clusters, of which, 14 were clearly ex- 
tended compared to the ROSAT point-spread func- 
tion. Combining this with results from H92, 40% 
of the z > 0.14 EMSS clusters are extended in the 
X-rays thus securing their classifications as X-ray 
clusters. For the other 7 clusters, only 2 were def- 
initely classed as point-like X-ray sources and the 
optical follow-up work of these clusters supports the 
idea that these maybe mis-classifications. The re- 
maining 5 clusters are uncertain, either due to a lack 
of signal-to-noise to make a conclusive statement on 
their extent, or, they simply resist easy categorisation 
(MS2137.3-2353, MS1512.4+3647). 

This exercise also highlights the power of the SHARC 
analysis methods. We were able to extract known 
distant X-ray clusters from the ROSAT PSPC data 
archive based upon their observed X-ray extent (clus- 
ters make up a small fraction of all X-ray emitting 
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sources). In cases where we did not detect an extent, 
there is a clear and well-understood reason for miss- 
ing the cluster. For example, certain systems (~ 10% 
of the EMSS) have intrinsically compact cores com- 
pared to typical X-ray clusters and the rest of the 
EMSS sample and may represent as a new class of X- 
ray source (see Stocke et al. 1991). Full details of the 
SHARC survey selection function will be presented in 
future papers. 

We have used this improved EMSS catalogue to re- 
visit the question of X-ray cluster luminosity evolu- 
tion. The original work of H92 suggested a steepening 
of the slope of the XCLF with redshift with a statis- 
tical significance of 3ct. We employed similar statisti- 
cal analysis as H92 (binned 1/V a test and Maximum 
Likelihood analysis) to our data and found that, in 
general, we find similar XCLFs as presented by H92. 
However, we found two significant differences. First, 
Figure 9 shows that the EMSS provides a poor deter- 
mination of the local X-ray cluster luminosity func- 
tion and shows that the low redshift EMSS XCLF is 
consistent with all the other EMSS XCLFs. Secondly, 
we find a shallower slope for the high redshift EMSS 
XCLF, which immediately reduces the significance of 
any claimed X-ray cluster luminosity evolution. 

We have compared our determinations of the EMSS 
XCLFs with those found in the literature. For this 
comparison, we have combined the two lower red- 
shift EMSS XCLF since there is little statistical evi- 
dence that these functions are significantly different. 
This combined low redshift XCLF (z = 0.14 to 0.3) 
agrees well with that derived from the RASS (Ebel- 
ing et al. 1995, 1996) and strongly suggests there is 
no evidence for XCLF evolution below z — 0.3. Fur- 
thermore, a 2D KS test between this combined low z 
XCLF and the high z XCLF indicates that the signif- 
icance of any evolution between the two is only 85%. 
If we compare the original low z EMSS shell (0.14 to 
0.2) with the high z shell, the significance of any ob- 
served evolution only rises to 89%. These findin gs are 
now consistent with those of Collins et al. 1997 . 

At greater redshift (z > 0.6), the data remains 
scarce, but at low cluster luminosities the XCLF does 
not appear to be radically different. At high luminosi- 
ties (£3,(0.3 - 3.5) > lO^ergs" 1 ), the fact that the 
EMSS does not detect such clusters may be indica- 
tive of evolution. Overall, Figure 10 summaries our 
new determinations of the EMSS XCLFs within the 
framework of other work in this field and, we believe, 
demonstrates that the significance of any claimed evo- 



lution with redshift can only be mild at best. Several 
projects are already underway (i.e. the SHARC sur- 
vey; Burke et al. 1995) which will help clarify the sit- 
uation by constructing larger samples of X-ray clus- 
ters from the ROSAT data (both RASS and archival 
data). 
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